Mon. Not. R. Astron. Soc. OOP. [TlfTOl COOOO') Printed 5 August 2008 (MN 1*1^ style file v2.2) 



Ellipsoidal collapse and the redshift space probability 
distribution function of dark matter 



oo 
O 

o 

<: 

in 



o 



(N 
> 
00 

m 
(N 

in 
o 
oo 
O 



Tsz Yan Lam* & Ravi K. Sheth* 

Department of Physics & Astronomy, University of Pennsylvania, 209 S. 33rd Street, Philadelphia, PA 19104, USA 



5 August 2008 



ABSTRACT 

We use the physics of ellipsoidal collapse to model the probability distribution func- 
tion of the smoothed dark matter density field in real and redshift space. We provide a 
simple approximation to the exact collapse model which shows clearly how the evolu- 
tion can be thought of as a mo dification of the s pherical evolution model as well as of 
the Zeldovich Approximation (IZerDovichlll97Clt l. In essence, our model specifies how 
the initial smoothed ovcrdcnsity and shear fields can be used to determine the shape 
and size of the region at later times. We use our parametrization to extend previous 
work on the real-space PDF so that it predicts the redshift space PDF as well. Our 
results are in good agreement with measurements of the PDF in simulations of clus- 
tering from Gaussian initial conditions down to scales on which the rms fluctuation 
is slightly greater than unity. We also show how the highly non-Gaussian non-linear 
redshifted density field in a numerical simulation can be transformed so that it pro- 
vides an estimate of the shape of the initial real-space PDF. When applied to our 
simulations, our method recovers the initial Gaussian PDF, provided the variance in 
the nonlinear smoothed field is less than 4. 

Key vifords: methods: analytical - dark matter - large scale structure of the universe 



1 INTRODUCTION 

The probability distribution function (hereafter PDF; some 
authors prefer to call this the probability density function) 
specifies the probability that a randomly placed cell of spec- 
ified size and shape contains a certain specified density. 
As gravitational instability alters the spatial distribution 
of dark matter, the PDF of the dark matter density field 
evolves. Broadly speaking, there are two approaches to mod- 
eling this evolution: one is to study the evolution of the lower 
order moments of the PDF, and the other studies the shape 
of the PDF in its entirety. 

Typically, the first approach attempts to describe the 
evolution exactly, using methods derived from perturbation 
theory, but then the PDF is often computed approximately, 
upon truncati ng higher o rder t erms in the perturba- 
tion theory ("Bcrnardeau' 'l994j; IScoccimarro fc FriemanI 
1 19991 : l^rnardcau ct al. 2002). The second uses sim- 
ple approximations to the exact evolution to derive 
tractable results; in almost all cases, these specify lo- 
cal transformations which describe how quantities in 
the initial field determine the evolved density. The best 
studied of thes e approxim a tions are the spherical evo- 
lution model l|Bernardeaul 1 19941 : IProtogeros fc Scherreij 
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19971 : iFosalba fc Gaztafiagal Il998l: iGaztanaga fc CroftI 
I999I : iBetancort-Riio fc Lopez-Corredoiral |2002|) . 



the Zeldovich approximatio n (iBetancort-Riic Il99ll: 



Padmanabhan fc SubramanianI Il993l : Kofman et al.l 



1994) and i ts extension to the ellipsoidal collapse model 



( Ohta et al] [2004: Lam fc Shoth 2008). In principle, there 
is a limit to how well ap proaches ba, sed on deterministic 
transformations can fare l|shethlll998l ). but, provided one 
restric ts attention to larg e enough scales, this is a small 
effect (|Lam fc Shethll2008l ). 

The associated PDF in redshift, rather than real space, 
is less well-studied. In this case too, the PDF is described in 
terms of its moments (|Hivon et al.lll99 F: 'Bernardeau et al 



2^002), using the spherical model (jScherrer fc Gaztana^ 



i 



2001), or using th e Zeldovich app r oxim ation (|Hui et al 
2000). In addition, IWatts fc Ta"^ (|200ll ) describe a cat 
culation of the PDF based on second order perturbation 
theory. The main goal of the present work is to present a 
calculation of the PDF which is based on the ellipsoidal col- 
lapse model. 

Our work is motivated by the fact that the spherical 
collapse model is inconsistent with a linear theory analy- 
sis of the second moment of the PDF by iKaiserl (|l987l ). 
Whereas Kaiser showed that the ratio of the redshift space 
variance to that in real space should be {1 + 2f /3 + /^/5) 
(here / = dln_D/dlna with D the linear theory growth fac- 
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tor and a the expansion factor), a finding which was con- 
firmed by iFisheif (|l995l ) using a very difi'erent approach, 
the spherical model predicts that this factor is (1 + //3)^ 
IScherrer fc Gaz tahaga 2001,) . An analysis based on the Zel- 
dovic h approximation does yield Kaiser's factor l|Qhta et alj 
I2OO4I ). providing yet another, very different derivation of 
this factor, so it is the spherical model which is in error. 
For this reason, our analysis of the redshift space PDF is 
based on the ellipso idal collapse model. Although this model 
has a long history llLin et al .11 19651 : Ilck3l 19731 : IWhite fc Sil^ 



19791: iBarrow fc Silk 



198l|), we use it in the form given by 



Bond fc Mverj ( 199€ ) . They showed how to write this model 



so that, to lowest order it reduces to the Zeldovich ap- 
proximation. It is this formulation we use, so the approach 
we outline below is guaranteed to correctly reproduce the 
Kaiser factor. This ellipsoidal collapse model (and its as- 
sociated tidal shear effects) has previously been used to 
motivate more accurate m odels of the abundance of viri- 
alized dark matter halos (jSheth et al.l I2OO1I ): it also pro- 
vides a framework for studying how the morphology, rather 
than just the density, of the la rge scale enviromcnt af- 
fect structure formation ( Shen et al. 2006; Dcsiacqucs 2008; 
iDesiacaues fc Smithlliooal v 

We describe our model in Section [S] and compare its 
predictions with measurements in a numerical simulation in 
Section [S] Section |4] describes a method for reconstructing 
the shape of the initial real-space PDF from the evolved 
redshift space PDF. A final section summarizes. Details of 
some of our methods are provided in two Appendices. 



2 THE REDSHIFT SPACE PDF 

In what follows, the variance of the initial density fiuctuation 
field when smoothed on scale Rm plays an important role. 
It is denoted by 

dk k^Pi^ik) 



aliM) 



27r2 



\W{kRi, 



(1) 



where Phik) is the power spectrum of the initial field, ex- 
trapolated using linear theory to the present time, W{x) is 
the Fourier transform of the smoothing window function, 
and Rm = (SM/Anp)^^^ . In addition, we use D{t) to denote 
the linear growth factor, and / = dlnD/dlna « 57^^®. 

2.1 The non-linear density in real space 

The nonlinear overdensity of a region of volume V containing 
mass M is 



1+5^—, 
pV 



(2) 



where p is the mean density. If the evolved Eulerian region 
is a sphere, then the evolved overdensity in real space is 
well-approximated by 



(1-^^/3)^ A 1 
= J-JTTs: 11(1 " A,) 



(3) 



(see Appendix . Here 5c is the critical value of spherical 
collapse model and its exact value depends weakly on the 
background cosmology. In a ACDM cosmology 5c ~ 1.66, 
while in an Einstein-de Sitter model 5c 



are the eigenvalues of the deformation tensor, extrapolated 
using linear theory to the present time, and 5i = Aj is 

the linear theory overdensity. 

This approximation has considerable intuitive appeal: 
the factor in the product sign is the Zeldovich approximation 
to the evolution, and the factor in front is the ratio of the 
evolution of a sphere in the Zeldovich approximation to that 
in the spherical collapse model. The expression actually uses 
a simple approximation to the evolution of the density in 
the spherical model; in principle, one could use the exact 
evolution instead. 



2.2 The non-linear density in redshift space 

The redshift space overdensity pn is related to the real space 
one by the mapping from real space to redshift space 



V ■ e ^ 
s = x~\ — e, 



(4) 



where s is the redshift space coordinate, x is the real space 
coordinate, v is the peculiar velocity at x, H is the Hubble 
parameter, and e is the line-of-sight direction unit vector. 
Then, 



ps — pr 



1 + 1^ 

H dxi 



(5) 



where we assume the line-of-sight direction is along the third 
axis. The velocity of the ellipsoidal region along its principal 
axis can be computed using equation (|A10[I . 



HRk 



= 1- 



/ [rI\, - A\5i [1 - (1 - 5i/&oY^"'-'] /3} 



Rid - A,) - [1 - 5,/3 - (1 - 5i/5cY^'^] ' 

(6) 

where / = dlnD/dlna and D{t) is the linear growth 



to A by equation (|A9|I . For a randomly oriented ellip- 
soid we label the Euler Angle components (61,62,63) = 
(cos i/) sin 0, sin 1/) sin 0, cos 0) for rotation from the line-of- 
sight coordinate to the coordinate of the principal axes of 
the ellipsoid. As a result the redshift space overdensity is 



(l-5;/3) 



{l~5i/5c) 



1-E 



/ {i?I,Afe - [1 - (1 - 5i/5cY^/^''-] 73} 



^ Rl{l - Afe) - Al [1 - 5,/3 - (1 - 5i/5cY^/^] 



6fe 



(7) 



Notice that pa pr when / 0, as it should. 



2.3 The PDF from ellipsoidal collapse 

Our redshift space perturbation theory calculation makes 
use of the same assumptions as in real space: there is a local 
mapping from the eigenvalues Ai of the deformation tensor 
to the non- linear overdensity p^; and the smoothing scale 
associated with ps in the initial field is the one containing 
the same mass with the correct shape. The most impor- 
tant difference in redshift space calculation is that the map- 
ping from \i to Ps is no longer deterministic: the orientation 
of the ellipsoid with respect to the line of sight introduces 
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stochasticity which must be integrated over. As a result, the 
evolved (Eulerian) PDF at fixed V is related to the initial 
(Lagrangian) PDF at fixed mass scale M by 

psP{ps\V)dps = y dAdep(A|cr)(5D [ps = pl(A,e)] , (8) 

where ps = M/M and p's(X,e) is the redshift space over- 
density given in equation ([7]). 

In p ractice the joint dis tribution oi p{\\a) is a function 
of Xi/a ()Doroshkevichlll97d l. where 



r 



k<j 



ln(# 

Rl 



(9) 



Here Rj den otes the initial size of the ,?'th principal axis , and 
B = 0.0486 ('Betancor t-Riio fc L6pez-Corredoirall2002l 'l. As 
a result, it is useful to rewrite equation ([8| as 



J dAdep(A|l) (5d [pa = \p's{(rX,e)\ 



Ps p{ps\V) dps 

m 

As happens for the real space calculation llKofman et all 



1 19941 : IScherrer fc Gaztanagall200l] : lOhta et al.i2004h . this es- 
timate of the redshift space PDF is not correctly normalized: 



N ; 



: / dpsP{ps) / 1. 

Jo 



(11) 



(On the other hand, the integral of equation [8] over all ps 
does equal unity.) Therefore, we set 

p's=Nps, and P? P {p's) = plp{Ps)- (12) 

This e nsures that f dp[, p'„ p' (p'„) and J dp'gp'{p'g) both equal 
unity (|Lam fc ShethluOOSi ). In practice, for the scales of in- 
terest in what follows, (cells of radius 4, 8, and 16/i~^Mpc), 
is only slightly larger than unity: 1.264, 1.104 and 1.052. 

2.4 Other approximations 

We will compare our model for the PDF with two other 
useful approximations which are less physically-motivated. 
One is the Lognormal, 



PsPips 



exp[-(lnps + psf/2a^] 



(13) 



where p,a = ul/^. and = In(p^ ). and cr^ is g iven by mul- 
tiplying the Kaiser factor by the lSmith et al.l (|2003l ) fitting 
formula for the real space variance a^. 

The second is based on the finding that, in simula- 
tions, the d istribution of {5a I (J a) is almost the same as that 
of {5r/(Jr) jScherrer fc Gaztafiagal l200ll ). Since the spher- 
ical model provides a good d escription of the [Sr/ur) in 
real space (|Lam fc Shethll2008h . we use it together with the 
Kaiser factor to set 



P{ps) = — Psc(pr), 
0"3 



(14) 



where Psc(p) is the real space PDF with the spherical col- 
lapse model. In essence, this can be thought of as a model 
in which the redshift-space overdensity is a linearly biased 
version of the real-space overdensity, the co nstant o f pro - 
portionality (the bias factor) being ("(7.,/o-r-). iKaiseij (|l987l ) 
suggests that this sho uld be an exc e llent approximation on 
large scales; although IScoccimarrol (|2004l ) suggests that it 



is not a well-motivated approximation on intermediate or 
small scales, measurements of the redshift space PDF by 
IScherrer fc Gaztafi agal (|200ll ') suggest that it is nevertheless 
reasonably accurate. 



3 COMPARISON WITH SIMULATION 

We compare our models with measur ements in a simula- 
tion which was run and analyzed by ISmith et al.l l|2007l ). 
The simulation box is a periodic cube 512/i~^Mpc on 
a side in a cosmological model where (fim, Oa, erg) = 
(0.27,0.73,0.72,0.9). The simulation followed the evolution 
of 400'^ particles, each of mass 1.57 X lO"/i-^M0. Halos 
were identified using the FOF algorithm with link length 
0.2 times the mean interparticle separation. 

To measure the Eulerian real-space PDFs we place 
spheres of volume V randomly in the simulation box, and 
record the number of particles in each sphere, taking care to 
account for periodic boundary conditions. We measure the 
redshift-space PDFs similarly, after transforming the parti- 
cle positions as follows. We assume that one of the sides of 
the box lies parallel to the line of sight. The redshift space 
coordinate is then given by adding the peculiar velocity (di- 
vided by the Hubble constant) in this direction to the real- 
space coordinate, taking care to account for the periodicity 
of the box. 

Technically, our model of the PDF is for a continuous 
density field rather than for discrete particle counts. In prac- 
tice, we are interested in large enough volumes that discrete- 
ness effects are negligible, so we simply set p= 1 + 5 = N/N . 
However, there is another effect which is more pernicious. 
The 'virial' motions of particles within halos give rise to 
long fingers-of-god — this makes halos stretch about 7 times 
longer al ong the l ine of sight than across — an effect first 
noted bv lJacksonl(|l972l ). which is entirely absent in the per- 
turbation theory calculation of the redshift space PDF. To 
illustrate their effect on the PDF, we have also computed 
redshift space positions after first setting all virial motions 
to zero, so that dark matter particles have the same speed 
as their associated halos. (I.e. we assign all particles in a 
halo the same velocity as that of the halo center of mass.) 
In the figures which follow, filled black squares show the true 
redshift space measurements, and open orange squares show 
results when fingers-of-god are absent. 

Figure [T] shows the redshift space PDF of the dark mat- 
ter in spheres of radius 8/i~^Mpc. The main panel shows 
this using a log-In scale to emphasize the behaviour of the 
low probability tails, and the inset shows the same data us- 
ing linear axes to highlight the behaviour near the peak of 
the distribution. Filled and open squares show the measure- 
ments in the simulation (i.e., before and after compressing 
the fingers-of-god). These indicate that fingers-of-god are 
most important in the denser cells; this is not unexpected, 
since dense cells host the most massive halos which have the 
largest virial motions. 

Our perturbation theory based prediction is represented 
by the dotted (blue) curve — notice that it provides a sig- 
nificantly better fit to the simulations when the fingers-of- 
god have been compressed. The solid (magenta) curve is the 
perturbation theory based prediction with i? = in equa- 
tion which corresponds to using our ellipsoidal collapse 
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Figure 1. Comparions of the measured PDFs with various models in 8/i~^Mpc sphere. The outer panel shows the log-In plot and the 
inner panel shows the linear plot. Filled (black) squares are measurements from the simulation, open (orange) squares are measurements 
with finger-of-god removed. Dotted (blue) curve is the prediction of the perturbation theory by applying equation ((7j. Solid (magenta) 
curve is the prediction by setting B = in equation JSJl. Dot-dashed (green) curve and long-dashed (red) curve are the lognormal and 
real-redshift space mapping empirial models respectively. 



model but ignoring the effect on the variance associated 
with differences from a spherical shape. The dotted curve 
{B = 0.0486) is always closer to the measurements than is 
the solid one (B = 0), cons istent with the real-space results 
presented in lLam fc ShethI ((2008 )), although the difference 
between the two curves is not large. In what follows, we 
only show results for the case when B = 0.0486. The dot- 
dashed (green) curve shows the Lognormal (equation I13|) 
and the long-dashed (red) curve shows equation (|14|l . Our 
model provides a better fit than these other more empirical 
curves, except in the highest density tail. 

Figure [2] shows similar results for spheres of radius 
4/i~^Mpc and 16/i~^Mpc. In the larger cells, fingers-of-god 



are irrelevant, and the Lognormal provides a slightly better 
description than does our model, which again overshoots at 
large densities. The approach which assumes a simple map- 
ping from real to redshift space overdensity fairs reasonably 
well at large densities, but underestimates the height of the 
peak of the distribution, and overestimates the PDF at lower 
densities. Our model is substantially better than the other 
two in the smaller 4/i~^Mpc cells. 
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Figure 2. Same as previous figure [T] but in spheres of radius Ah ^Mpc (left) and 16/i ^Mpc (right). The prediction with spherical 
variance (the solid (magenta) curve in figure [ij is not shown. 



RECONSTRUCTION OF THE INITIAL 
DISTRIBUTIONS 



In|L am fc ShethI \ 20081 ) we described how the spherical col- 
lapse model could be used as the basis of a method for re- 
constructing the shape of the initial PDF from the Eulerian 
one. In principle, the method can be used as a consistency 
check for Gaussian initial conditions. We showed that, in 
real space, it can reconstruct the Gaussian form of the ini- 
tial distribution from the highly skewed nonlinear PDF on 
scales where the rms fluctuation is about 2 or less. 

Extending this reconstruction method to redshift space 
is complicated by the fact that we know the spherical model 
is inadequate (it does not reproduce the Kaiser factor). How- 
ever, the ellipsoidal collapse model, which we have shown 
works quite well, has many more parameters — the lengths 
and orientations of the three principal axes — so that the 
final density is determined by a stochastic rather than de- 
terministic mapping. 

So we have tried a simpler approach which is based on 
the fact that equation (|14|l provides a reasonable description 
of the redshift space PDF. Namely, we assume that 

5s = 5r {(Jslfyv), (15) 

and use the (square root of the) Kaiser factor in place of 
(Js/<jr', we have found that using our ellipsoidal collapse cal- 
culation of this quantity instead produces slightly worse re- 
sults. Then we set 

l-[(l-f5.a./a.)/jV.c]-^/^° ^ (^g^ 



<7L(A//iV,c)/<5c 

where 5c is the critical linear density associated with collapse 
in the spherical model, and N^c is the normalization factor 
(the analogue of equation 1 in the spherical model calcu- 
lation of the real space PDF. Finally, we make a histogram 
of u by weighting each cell which contributes to p{5s) by its 
value of (1 -I- Ssar/as)/Nsc- 

Before showing how well this works. Figure |3] shows a 
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Figure 3. Comparsion of real- and redshift-space overdensities in 
cells of radius 4/i~lMpc (top), 8/i~lMpc (middle) and 16/i~lMpc 
(bottom). Filled and open symbols are before and after fingers of 
god have been removed. They show the mean, and error bars show 
the rms, for narrow bins in {S/a)^. If equation I I15I I is accurate, 
then {S/a)s — {&/a)r = with little scatter. 



test of our assumption that equation (I15|l is accurate, mean- 
ing that 5s is linearly related to 5r , with little scatter around 
the mean relation. Evidently, these assumptions are reason- 
able on large scales (bottom and middle panels) , but become 
increasingly worse on smaller scales (top). 

Figure |4] shows the reconstruction results for 4, 8 and 
16/i~^Mpc spheres. The left- and right-hand panels show the 
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Figure 4. Comparsion of the reconstructed linear PDF of <5/ct 
(histogram) with the expected zero-mean unit-rms Gaussian form 
(solid curve). Upper panels show results of reconstruction for 
4/i~^Mpc scale, middle panels show results for 8/i~^Mpc. lower 
panels show results for 16/i~^Mpc. Left- and right-hand panels 
show results of reconstructions which begin from the filled and 
open squares in the previous Figures (i.e. estimates oipi&s) with 
and without fingers-of-god) . Dotted curves in each panel show the 
corresponding nonlinear PDFs (equation (S)! . 

result of transforming the PDFs shown by filled and open 
squares in the previous Figures (i.e., full and compressed fin- 
gers of god) . Smooth solid curves show a Gaussian distribu- 
tion for comparison. The reconstruction method works well 
for the 8 and 16/i~^Mpc scales. This method even works for 
the smaller 4/i~^Mpc scales if the fingers-of-god have been 
removed. The dotted curves show the corresponding non- 
linear distributions (equation IH} to demonstrate that our 
reconstruction method transforms significantly skewed dis- 
tributions back to the original symmetric gaussian shape. 



5 DISCUSSION 

We have extended our previous work on the real space dark 
matter PDF to redshift space. To do so, we provided a sim- 
ple parametrization of the evolution of the nonlinear density 
in the ellipsoidal collapse model. This shows explicitly how 
the model extends both the spherical collapse model, and 
the Zeldovich approximation (equations lAlOl and (Tjl . This 
parametrization, with our prescription for relating spatial 
statistics in the initial and evolved fields (equation[8]), results 
in excellent agreement with measurements from a numerical 
simulation, in both real (Figure FAlll and redshift space (Fig- 
ures [T] and [Sjl. Our parametrization fares significantly better 
than analyses which are based on exact second and third 
order perturbative expansions of the ellipsoidal collapse evo- 
lution — these produce PDFs with a power-law tail at high 
densities which is a signature that the analysis is breaking 



down — a similar tail was seen for the Zeldovich approxima- 
tion (|Hui et al.ll2OO0l ). Our approach also fares better than 
the Lognormal model, and one which assumes that redshift 
space overdensities are linearly proportional to those in real- 
space: 5s — {(Js/or) 5r where the constant of proportionality 
is given by the (square-root of the) ratio of the variances. 
The latter is known to be a good approximation on large 
scales — we find that it is accurate on small scales, except 
in underdense regions. 

We also used our results to motivate a method for re- 
constructing the shape of the initial PDF from the nonlin- 
ear evolved redshift space PDF. This method works well on 
scales where the rms fluctuation is less than 2 when fingers- 
of-god have been removed (Figure |4]). We are in the pro- 
cess of applying this method to reconstruct the baryonic 
ascoustic oscillation (BAO) peak in the correlation function 
(Eiscnstcin ct al. 2007). Instead of using equation (I16|) it 
is also possible to formulate the reconstruction using equa- 
tions (|A10|) or (O; this is also the subject of work in progress. 

While this provides a nice graphic reconstruction of the 
Gaussianity of the initial conditions, more work is neces- 
sary before it can be used to distinguish between a truly 
Gaussian model and the very mildly non-G aussian models 
which are not yet excluded by the CMB (|Komatsu et al.l 
2008). In this respect, it may be interesting to extend our 
method for estimating the nonlinear PDF to the family of 
non- Gaussian models which are currently studied in the con- 
text of using cluster abunda nces and clustering to constrain 
primordial non-Gaussianit y jSlosar et al.l l2008l : iDalal et al.l 
l2008l ; iLo Verde et"ai]|2008l ) . 
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APPENDIX A: RELATION BETWEEN INITIAL 
AND FINAL DENSITY IN THE ELLIPSOIDAL 
COLLAPSE MODEL 

In this Appendix we derive relatively simple analytic ap- 
proximations to the evolution of the density within a triax- 



ial perturbation. The exact evolution of the perturbation is 
as sumed to be describe d by the ellipsoidal collapse model 
of lBond fc MversI l| 19961 ). IShen et al.( l|2006l ) provide a sim- 
ple but useful approximation to this evolution which is based 
on a beautiful analy sis of the collapsing ellipsoid problem by 
IWhite fc Sill3 l| 19791 ). It is this approximate solution which 
we exploit. 



Al The collapse model 

Gravity evolves the initial axes {R\, R\, R\) of a triaxial 
object to (Ri, R2, R-j,) because 



1 + S 



— -t Aext.fc 



(Al) 



dt2 u - " \^ 3 2 

where the interior and exterior tidal forces are respectively, 

dr 



/2 



-^oxt,fc(^) — 



2 

^ 3' 

m 

D(U) 



(A2) 
(A3) 



teond fc Mverslll996l ). Here Afc are the initial eigenvalues of 
the deformation tensor, and the initial conditions of equa- 
tion (|A1|) are set by the Zeldovich Approximation: 



Rk{ti) ^ a,[l ~ \k{U)] 
dRk{U 



dt 



H{U)Rk{ti) - a{U)H{U)fi\{ti), 



(A4) 
(A5) 



where /i = d In D/d In a and D{t) is the linear growth factor. 

Altho ugh equation (lAll) in general must be solved nu- 
merically, IShen et al.l l|2006l ) showed that the following ap- 
proximation was rather accurate: 



Rk(t)=^Rk{U)[l-D{t)\k\ 



ait) 



Rh{ti) 



Dit) a,{t) 



-5r 



a{t) 



(A6) 



where Rh{ti) — i/ Rj{ti)~^ and ae{t) is the expan- 
sion factor of a spherical universe with initial overdensity 
Si = ^j{ti)- This ex pression is inspired b y an analysis of 
collapsing elhpsoids bv lWhite fc Silkl l| 19791 ). The first term 
in the expression above is the Zeldovich Approximation to 
the evolution; the second term is a correction to this approx- 
imation. When A2 = Si/ 3 then the evolution of the second 
axis is almost exactly like that of a sphere with initial over- 
density Si. 



A2 Evolution to a spherical volume 

For an Eulerian sphere of Re, equation HA6|) implies that 

R[[l - D{t)\^] = Rl[l ^ D{t)\2] ^ Rill - D{t)\3]. (A7) 
Substitute into the definition of Rh{ti) 

3 _ Rk(U)(l- D{t)Xk) 



Rh{U) = 



1 - D{t)Sr/3 



(A8) 
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A3 Evolution to an arbitrary shape 




Figure Al. The real space non-linear PDF. Symbols show mea- 
surements in a numerical simulation (from Lam & Sheth 2008). 
The solid (blue) curve uses our approximation to ellipsoidal col- 
lapse (equation I AlOl l and the dashed (red) curve is based on our 
approximation to spherical collapse (equation IA11[| . 



where = 1, 2, 3. Therefore we can write the evolved Rk{t) 
as 

Re _ ae{t)/a(t) 
Rl' 1- Dm/3 



[1 - I)(i)Afc] 



(A9) 



where Kl = Rkiti). This writes the evolution as that in the 
Zeldovich Approximation, times a correction factor which is 
simply the ratio of the evolved and initial radii of a Zeldovich 
sphere (1 — DSi/3) divided by the same ratio in the exact 
spherical evolution model. Thus, the evolved overdensity is 



1 + 5 



3 ni 



(1 - D{t)S,/3) 



Re 

{1-D{t)5,/3Y 



,=1-^ KW/a(t)) 

3 3 ^ 

(l-D(i)5./5.)^ nT~D(t)A, 



3 



(AlO) 



The final expression uses a simple approximation for the 
evolution in the spherical collapse model: 



ae{t) 



Oc 



-«c/3 



(All) 



iBernardeaull 19941 : IShethlll998l ) . where 5^ « 5/3. Notice that 
when all the eigenvalues are identical, then equation (|A10|) 
reduces to the spherical collapse model. 

Insertion of equation (lAllfl or (IA10[) into equation ([8]) 
of the main text yields models for the redshift-space PDF as- 
sociated with spherical and ellipsoidal collapse, respectively. 
The corresponding models for the real-space PDF are given 
by setting / = in that expression. These real-space mod- 
els are shown by the dashed (red) or solid (blue) curves; 
both are in excellent agreement with the measured real- 
space PDF (symbols). 



When the Eulerian volume in consideration is triaxial, the 
computation of the overdensity is less trivial. In general 
one must solve a system of non-linear equations for = 
a{t)Rl/a,Ru{t) in 



Xk[l-D{t)Xk\ = 



3 1 



3 a{t) 



E 



-Rfe_l_ 

Ri xi 



(A12) 



The overdensity is l-|-(5 = 11^= i ^j- The above equation can 
be solved numerically or iteratively with 



JO) 



1+3 1- 



3 a{t) 



E Rk 1 



[1 - D{t)\k] 



a,{t)/a{t) l-D{t)\k' 



(A13) 
(A14) 



An analytical expression for the overdensity of the ellipsoid 
can be obtained if we make the approximation that the col- 
lapse o f the second axis i s very close to that in the spherical 
model. IShen et all l|2006h show why this is almost always an 
excellent approximation. One can then approximate the col- 
lapse of the second axis using either of the two expressions 
below: 



Rh 



l-D{t)S,/3 



1 



R2{t) 



D(t)5, 



5c 



-«c/3 



(A15) 



The first expression is our approximate description of evo- 
lution to an Eulerian sphere, and the second is simply the 
spherical collapse model. 

Substituting these approximations into equation (|A6P 
yields 



||[l-i3(t)A,]= <i 



Rk _ ^ 
Rk 



R2 



- 1 



1 - g(t)g,/3 
ae{t)/a{t) ' 

l-D{t)\2 

aeit)/a(t) ' 



(A16) 



for which the associated overdensities are 

3 

P ^ 



= < 



R, 



n 

j=i 

3 

n 



a{t) l-D{t)Si/3 
ae{t) l~D{t)\j 

a{t) 1-D{t)\2 



R2 a^it) 
R, 



R2 1 Rj 



Ra 
R, 



a{t) l-D{t)Si/3 

ae{t) 1 - R2/RJ 
a{t) 1-D{t)\2 



(A17) 



where ae(t) is approximated by equation (|A11|I . Compared 
to equation (|A10[) the last term in the first approximation 
is the extra factor accounting for the non-spherical shape in 
the final volume. For Ri = R2 = R3, this extra factor is 
unity and equation (|A17P has the same expression as in the 
spherical case. 



Non-linear dark matter PDF in redshift space 9 



APPENDIX B: PERTURBATIVE TREATMENT 
OF REDSHIFT SPACE DISTORTIONS 



iKaiseil (| 19871 ) derives a simpl e relation betw een the real and 
redshift space power spectra. iFisheil (|l995l ) provides a very 
different derivation of Kaiser's expression. The Zeldovich ap- 
proximation prov ides yet another derivation of this relation 
l|Ohta et al.ll2004l '). We show that our approximation for el- 
lipsoidal collapse also reproduces this expression, and then 
use it to deriv e the next or der corrections. This analysis is 
useful because IScoccimarrol (2004) has shown that Kaiser's 
approximation is accurate only on rather large scales. 

Before we present the algebra, it is worth noting that 
the approximation used in the main text for the real space 
evolution of a collapsing ellipsoid follows from (i) writing 
the full ellipsoidal collapse as the Zeldovich approximation 
times a correction factor which is the ratio of the evolution 
of a sphere in the Zeldovich approximation to that in the 
spherical collapse model; and (ii) using our convenient ap- 
proximation to the exact evolution in the spherical collapse 
model. In the perturbative analysis which follows, we have 
chosen notation which illustrates what happens if we keep 
approximation (i), but treat the spherical collapse model 
exactly rather than approximately. 

We begin by writing the relation between the linear and 
nonlinear overdensities in the spherical evolution model as 

a{t) - ^ 



Then the exact result has 



1 + E^^"- 



(Bl) 



1^2 



92 

9l 



10 



12^ + ^+6^ 

gi 91 9i 



(B2) 



where gi is the ith order growth factor (e.g. Fos- 
alba & Gaztaiiaga 1998). For A / (32/g?) ~ 
-(3/7)J7-i/i«, (.a3 c./g?) ^ -(l/3 )n-^/f^^ and {gz./gf) ^ 
(10/21) iBouchet et al] (Il995l ). The dependence 

on SI is so weak that we can use the S7 = 1 values to find 
1^2 ~ 1.62 and ~ 3.926. In contrast, our previous approx- 
imation has 



n 



(B3) 



so i/2 ~ 1.6 and ~ 3.52. 

With this notation, our approximation for the evolution 
of an ellipsoid is 

3 / 00 \ 3 



l + 5r= [1 



/ 00 \ 3 

\ n=l ' / k-1 



(B4) 



where 



^ (2) 1^2 .2 ,Sf . , 



2 ' 3 



^V*"' = ^ <5f + ^ 5? - 25, ^ A, Afc + Al A2 A3 . (B5) 



In each expression above, the first terms are those associ- 
ated with the spherical model, and the others are the cor- 
rections which come from the ellipsoidal collapse. The terms 
associated with a perturbation theory analysis of the exact 
ellipsoidal collapse model are given by Ohta et al. (2004). 
The expressions above show the terms which are associated 
with our simple approximation to the exact evolution. They 
differ slightly from those associated with the exact analysis; 
e.g., the second order expression associated with the exact 
analysis is (4/17)((5f — 3/2), but they otherwise have the 
same form, suggesting that our simple expression captures 
the essence of the ellipsoidal collapse evolution. 
The overdensity in redshift space is 

l+5r 



=l + 5«+5?'+5f' + 



(B6) 



where 



9k = 



>^ - [fiSi - (1 + Er=i i^sr)-'^' (E™.i u^sr)] 73} 



Rlil - Xk) - Al[l - 5,/3 - (1 + Er=i 



-1/31 



pi)Sf + fiXl 



+ 



,V2 4 !/3 „ 2U2 ^'^■2f,r,f \X 

(y ,/2 - 3.;i)Afc + {—fa ^fi ^12 + 2fi)di 



|/i)Afe5f + /iAfe 



fc -r 1 ^ g 9' 



(B7) 



and fi = din gi/ da. Thus, 



(B8) 



and 
Ai^' 

A (2) 
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k Bk 



2 /i 



5f 2 

-TT fife 



+ /i ^ Afc e| -I- /i ^ Aj Afe e]e\ 
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=.aE 



1^ f2_ 
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2z^ 
3 



5t^2 h 
6 /i 



+ 2) Si 



2 



1^2 f2 



2 /i 



2 2 

€,j Cf^ 



(B9) 



+ 2/f 

3 

+ /i ^ AjAjAfc e^e^efe. 

i,j,k—l 

Our leading order term, ^i^' is the same as that of Ohta 
et al. (2004), who showed that = {(Si^^)^), where the 
angle brackets denote the result of averaging over the Euler 
angles as well as the distribution of the Ai, leads to Kaiser's 
formula. The logic which leads to our second order term is 
consistent with that of Watts & Taylor (2001), except that 
because our collapse model considers motions relative to the 
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Figure Bl. Four approximations to the redsiiift space PDF. Solid 
line shows the prediction associated with the Zeldovich approxi- 
mation, but assuming that the associated cr is the same for all cells 
of volume V (Hui et al. 2000). Dashed and dot-dashed curves show 
the first- and third-order exact ellipsoidal collapse based analyses; 
these do account for the fact that cr depends on the Lagrangian 
size and shape (e.g. Ohta et al. 2004). Dotted curve shows the 
PDF associated with our approximation to the exact evolution. 

center of the object, it has no term which accounts for the 
motion of the center of mass. There are additional small 
differences which arise from the fact that our description 
of ellipsoidal collapse is based on our approximate model - 
had we computed our perturbation series expansion using 
the exact collapse model, we would have reproduced their 
expressions. 

Figure IBll illustrates one reason why we have used 
our approximation for the nonlinear dynamics, rather than 
worked with exact perturbative expansions to higher and 
higher orders. The dotted line shows the redshift space PDF 
associated with our approximation; the main text shows that 
it provides an good description of the simulations. The solid 
line shows the result of using the Zeldovich approximation 
for the evolution, and setting cr = ai pV) whatever the shape 
and overdensity of the initial object iHui et all ()2000l ). This 
results in a PDF with a tail that scales as pj^ at high den- 
sities — a signature that the approach has broken down. 
The dashed curve shows the result of using equation ((Oj for 
cr instead. This helps somewhat — although the onset of 
the pj"^ tail is shifted to larger ps, the shape of the PDF 
in underdense regions is adversely affected. The dot-dashed 
curve shows the result of going to third-order in the exact 
ellipsoidal collapse model (recall that, in this context, the 
Zeldovich approximation is like the first order term). This 
modifies the shape of the PDF slightly, but the problem at 
small ps remains. 



